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ENTROPY ANALYSIS OF KINETIC FLUX VECTOR SPLITTING SCHEMES 
FOR THE COMPRESSIBLE EULER EQUATIONS * 

SHIUHONG LUF AND KUN XU* 


Abstract. Flux Vector Splitting (FVS) scheme is one group of approximate Ricmann solvers for the 
compressible Euler equations. In this paper, the discretized entropy condition of the Kinetic Flux Vector 
Splitting (KFVS) scheme based on the gas-kinetic theory is proved. The proof of the entropy condition 
involves the entropy definition difference between the distinguishable and indistinguishable particles. 
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1. Introduction. There are many numerical approaches to the solution of the Euler equations. Go- 
dunov and Boltzmann schemes arc two of them [4]. Broadly speaking, Godunov scheme is based on the 
Riemann solution in the gas evolution stage, and the Boltzmann scheme uses the microscopic particle distri- 
bution function as the basis to construct the fluxes. While the construction methodology is different between 
the Godunov and kinetic schemes, both first order schemes can be written in the framework of the 3— point 
conservative methods. 

There arc mainly two kinds of gas- kinetic schemes, and the differences arc in the governing equations 
in the gas evolution stage. One of the well-known kinetic schemes is called KFVS which is based on the 
collisionless Boltzmann equation [9, 10], and the other is based on the collisional BGK model [15]. By 
combining the dynamical effects from the gas evolution stage and projection stage, the real governing equation 
for both KFVS and BGK schemes arc physically the same except the particle collision time r in the BGK 
scheme is replaced by the CFL time step At in the KFVS scheme [14]. 

The previous paper [ 11 ] analyzed the positivity property, such as positive density and pressure, for the 
gas- kinetic scheme. In this sequel, we analyze the entropy condition for the first order KFVS schemes. 

2. Preliminaries. We consider the one dimensional Euler equations of gas dynamics: 


(2.1) 


Pt +m x = 0, 
m t + (mU + p) x = 0, 
k Et + {EU 4- pU) x = 0 , 


where p is the density, U the velocity, m = pU the momentum, E = \pU 2 -f pe the energy per unit mass, e 
the internal energy density, p the pressure. We assume that the gas is a 7 -law gas, i.c., p = ( 7 — l)pe. In order 
to obtain the approximate solution for the above equations, the gas-kinetic scheme solves the Boltzmann 
equation in the gas evolution stage. 
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Council through RGC HKUST6166/97P and HKUST726/96E. 
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The Boltzmann equation in the 1-D case is [6] 


ft + ufx — Q(f, f), 

where / is the gas-distribution function, u the particle velocity, and Q{f, /) the collision term. The collision 
term is an integral function which accounts for the binary collisions. In most cases, the collision term can 
be simplified and the BGK model is the most successful one [1], 

Q(fJ) = (9-f)/r, 

where g is the equilibrium state and r the collision time. For the Euler equations, the equilibrium state g is 
a Maxwellian, 

JC + 1 

( 2 . 2 ) = 2 

where £ is a K dimensional vector which accounts for the internal degrees of freedom, such as molecular 
rotation and vibrations, and £ 2 = £? + £| + ■■■ + &■ Note that K is related to the specific heat ratio 7, 

K = (3 -7)/(7 - 1). 

Monotonic gas has 7 = 5/3, and diatomic gas has 7 = 1.4. The lower limit of 7 is 1, which corresponds 
to an infinite number of internal degrees of freedom. For example, 7 = 103/101 is equivalent to K = 100, 
which gives 98 internal degrees of freedom for the molecule. In the equilibrium state, A is related to the gas 
temperature T 
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where the entropy density is 


H= / /In fdudt 


and the corresponding entropy flux is 


G = 


j uf in fdud£. 


The first-order numerical conservative scheme can be written as 




■ * 7 + 1 / 2 ). 


where Wj ~ (pj^rrij, Ej) T is the cell averaged conservative quantities. F 7 4l/2 is the corresponding fluxes 
across the cell interface, and a = At/ Ax. For the lst-ordcr gas-kinetic scheme, the numerical fluxes across 
cell interface depend on the gas distribution function /j l + i /2 v * a (2-3). The discretized entropy condition for 
the above 3-point method is 


(2.4) 


H n+l < jjn 


+ ^ a - 


l/2-C" +1/2 ), 


where Hj = f f 3 In fjdud£ is the cell averaged entropy density and Gj+i/ 2 = / uf J+ i/ 2 In f 3 +\/ 2 dud£ > is the 
entropy flux across a cell interface. In this paper, wc prove the above inequality for the KFVS scheme. Since 
the KFVS scheme assumes an equilibrium distribution inside cell j at the beginning of each time step, H™ 
becomes 


= j <?" In g™ dudf; 

K -}- 1 

(2.5) = p™ In p n 3 + pj — — — (In — — 1). (with the equilibrium distribution in Eq.(2.2)) 

Since at the beginning of each time step, the gases in the cells j — 1, j, and j- f* 1 arc basically distinguishable, 
the updated flow variables WJ 1 ^ 1 inside cell j at time step n-f- 1 are composed of three distinguishable species 
from cells j — 1, j, and j + 1. So, the total entropy density is the addition of the entropy of different 

species. 

It is very difficult get a rigorous proof of the discretized entropy condition (2.4) for the nonlinear hy- 
perbolic system. The difficulty is mostly in the interaction between numerical gas from different cells. The 
update of the entropy in each cell is a complicated function of all flow variables including the ones from the 
surrounding cells. Since the entropy condition only tells us the possible direction for a system to evolve, it 
docs not point out exactly which way to go. So, in order to analyze the entropy condition for the discretized 
scheme, we design a “physical path” for the gas system to evolve. With the same initial and final conditions 
for the mass, momentum and energy inside each cell, the proof of the entropy condition becomes the proofs 
of the entropy-satisfying solution in each section of the physical path. Fortunately, for the KFVS scheme, we 
can design such a physical process. To show (2.4), we use results in statistical mechanics about the definition 
of entropy for distinguishable and indistinguishable particles. 


3. KFVS Scheme. In this section wc consider the kinetic flux-splitting scheme (i.e. collisionlcss 
scheme) proposed by Pullin [10] and Deshpande [2]. The scheme uses the fact that the Euler equations 
(2.1) arc the moments of the Boltzmann equation when the velocity repartition function is Maxwellian. As 
numerically analyzed in [7], the flux function of the KFVS scheme is almost identical to the FVS flux of 
van Leer [13]. In Section 3.1 we briefly recall the collisionlcss scheme. In Section 3.2 wc prove the entropy 
condition for KFVS under the standard CFL condition. The positivity of the KFVS scheme has been 
analyzed in [3, 9, 11]. 


3 



3.1. Numerical scheme. In order to derive the collisionlcss Boltzmann scheme, we need to construct 
the numerical fluxes across each cell interface. We suppose that the initial data (p(x),m(x), E(x)) arc 
piecewise constant over the cells Cj = [x j _ 1/2 ,x j+1/2 ]. At each time level, once Pj,m 3 and E } are given, the 
corresponding Uj and Xj can be obtained by the following formulae: 

(3.1) m = pU, E = \pU 2 + 

2 4A 

Let 

K ± 1 

(3.2) g . = p .^\ 2 e -x A( u-u 3 f+e) 

be a Maxwellian distribution in the cell Cj. The corresponding distribution function at the cell interface is 
defined by 


(3.3) 


f{ X i+l/2,t,U,Q 


{ gj, if u > 0 
9j+ 1, if w < 0. 


Using the formulae (2.3), we obtain the numerical fluxes 


(3.4) 


Fp.j+\/2 < 

F m,j+ 1/2 
F E,j+ 1/2 / 

( 


( ^crfcC-^Aj Uj)+¥^± 


= Pi 


(if + 4Jj) erfc(- y/XjUj) + 


j J_ 3 3 




+ Pj + 1 


— crfc( y/\j^~iUj + i ) - 




(% + 4 A^r) erfc( v A^7f/ >+1 ) - 


A i+ lU i+i 


v^7 


I K"+3 T T 




/ 


where the complementary error function, which is a special case of the incomplete gamma function, is defined 

by 


2 Z* 00 2 

erfebr) = — p= / e _t dt. 

V* Jx 

Using the above numerical fluxes, we are able to update p^m^Ej with the standard conservative formula- 
tions: 


(3.5) 



where W 3 = tU” +1 and 


( Pi ^ 

( Epj- 1/2 - Fp,j+ 1/2 \ 

rrij 

F 17 1 F m,j- 1/2 ~ F m j + i/2 1 

Ej y 1 

\ F E,j-\/2 ~ F E j+ i/ 2 ) 


A* 

^ ~~ Ax 5 

with A t the stepsize in time, and Ax the mesh size in space. The scheme can be viewed as consisting of the 
following three steps (although it is not typically implemented this way): 
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ALGORITHM (KFVS Approach) 

1. Given data {p™, C/j 1 , £”}, compute {A™} using (3.1). 

2. Compute the numerical flux {F p j+i/ 2 i ^m,j+i/2> ^Ej+1/2} using (3.4). 

3. Update using (3.5). This gives {p™ +1 , ra" +1 , }. 

3.2. Entropy analysis. The analysis of entropy condition for the KFVS scheme has attracted some 
attention in the past years. In [2], Deshpande stated the entropy condition in the smooth flow regions. In 
[5], Khobalatte and Pcrthame gave a proof of the maximum principle entropy condition for a gas kinetic 
scheme with a specific equilibrium distribution and a piecewise constant entropy function. In [8], an entropy 
inequality is introduced for a special distribution function. In this section, for the first time, we show that at 
the discretized level, the KFVS scheme satisfies the entropy condition with the exact equilibrium Maxwellian 
distribution. 

With the same initial and final mass, momentum and energy densities in Eq.(3.5), we can design a 
physical path for the flow updating process. The proof of the entropy condition is based on the entropy- 
satisfying solution in each section of the evolving path. 

In the first step, we consider the case when there is only gas flowing out from cell Cj. This gives 


(3.6) 


W* = 


Pj \ 

m* 

E 'i ) 


Pj 

= I rm 


+ <7 


/„<o ugjdudt - f u>0 ugjdud( 

/„< 0 u2 9 jdud£ - J u>0 u 2 gj dud{ 

V /„< 0 f (“ 2 + ¥) 9 idud£, - / u>0 § {u 2 + i 2 )gjdudZ 


The second step is to consider the inflow from adjacent cell Cj- 1, 

Iu> 0 m-idudt 
f U 2 Qj-\dud£ 

Lo%(v 2 + ?)9> -idudZ 




( 

Pj ^ 

\ ( 

(3.7) 

w = 


nij 

= a \ 



{ 

Ej ) 

1 \ 


In the third step, the inflow from adjacent cell Cj+\ is considered, 


(3.8) 



Pi > 


( ~ f u <o u 9j+idud£ \ 

W - 

rhj 

I = <7 

~ J u <o u2 9j+idudt 


Ej ) 

< 



The fourth step is to include particle collisions to let W* , W and W in the above equations to exchange 
momentum and energy inside cell j and to form the individual equilibrium states W * r , W f and W' with a 
common velocity and temperature, 




f Pi 1 


Pi \ 

Pi \ 

Pj \ 

Thj 


m i + 

rhj + 

rhj 

Ei J 


e: / 

Ei ) 

\ E; / 


(3.9) 


Pi 


m- 


V < / 



During the above collisional phase, the individual mass, total momentum and energy are unchanged. It 
can be verified that (pj,fhj, Ej) obtained by (3.9) arc exactly the same as those obtained by using (3.5). 
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In terms of updating conservative variables, the above four stages form the complete KFVS scheme. The 
entropy density H" + at time n + 1 inside cell C is the sum of the individual entropy of different species. 
Suppose that the CFL condition 


(3.10) 


a < 


max, {{Ujl + Cj) 


is satisfied, where Cj — 7 / 2 A„ is the local speed of sound. It has been shown in [11] that the positivity 

conditions arc precisely satisfied for the flow variables p* > 0 and pjEj - ~(mj) 2 > 0, as well as pj > 0 and 

pjEj - ^(rhj) 2 > 0. 

In the following, we prove that the discretized entropy condition is satisfied in the above four physical 
processes. As a result, the whole numerical path in the flow updating scheme satisfies the entropy condition 
(2.4). 

Lemma 3.1. Assume that the CFL condition is satisfied. If pj > 0 and PjEj > ^m 2 , then the entropy 
condition is satisfied in the updating process for [pj. 

Proof. Wc need to show that 


(3.11) f gj In gj dudS, < f g 3 Ingjdud^ + a f ugj In gjdud£ — f ugj In (/, 

J ~°o Jog lJu<0 Ju>0 

Wc use the following relations to express the * states in terms of the j states. 

^2^ jaj + j > 


dud(j 


Pj = Pj ~ a Pj 
mj = mj — apj 

E] = Ej - o Pj 


U 2 

2 + 4A 


IP 




A + 3„ \ 

+ IaT^ ^ + 


U]_ K + 2 
2 + 4A 


)4 


where 

(3.12) 


Q, 


erfc (-yfcUj) - erfc (i /XjUj') ; fy 


e~ x * u f 

\/*\j 


The equihbrium state gj has an Maxwellian distribution which corresponds to the macroscopic densities 
{pj,mj,Ej). 

After some algebra, 


I 9j lng*dud£ - / g 3 In gjdud£ - o / ug 3 In g 3 dud^ - / ug 3 In g } dud^ 

J -oo l-oo L/u<0 Ju> 0 


= PjF, 


where 
F - 


{ (l - \Wm+V>,)) ((A- + 2) In (l - + 2ft)) Inft) - fft}. 


a Xi 


K- 1-1 


*■ - 1 - + + ^) - (> - f + w) { {«! + *£) "«+ 


(2u; 


2 + A + 2 


Pj f + 


K +\ { ( U I + 2\j ) ajUj + 2U ^ j } 


2«t 2 A 
A + 


»{( 


Uf J_’ 
2 + 4Aj 


a j + Ujfij 
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The goal is to show that F < 0 for all positive a up to the CFL limit. We can reduce the number of 
parameters by one by introducing the non-dimensional number 2 = \/XjUj which is equivalent to the local 
Mach number in cell j. Wc also replace the parameter a by c € ( 0 , 1 ] (CFL number) which is defined by: 


( 3 . 13 ) 


C \f^3 

\A + \/ t 72 


Let 


Then 


OLiZ e * , c 

<t> = —— H 7 — = erf(z) z + 

1 y 7T 

e~ 2z2 >c 

ip — h z crf(z' 


A’ 

crf(z) 2 




v/t 72 ' 


F=(l-d<p) 


ln(l — d<j>) - ln/i 


2s/w 


where 


h= 1 


d 

(AT 4- 1)(1 -#) 2 




We now proceed to show that A = jF(z, A^, c) < 0, where the arguments of the function arc related to 
Mach number , gas constant, and CFL number, respectively. First note that F is an even function of z and 
hence we can restrict to the case z > 0 . By a direct calculation, d<j)/dz = erf(z) > 0 for 2 > 0 and thus <j> is 
minimum at 2 = 0 where it equals 1 /y/n. This shows that <f> is a positive function. 

Next wc show that 1 — is positive and less than one. Since both d and 4> arc positive, it is clearly less 
than one. To show that it is positive, it is sufficient to show this for c — 1 . Noting that c~ z /y/n < l /\/2 , 
we have 


0 < 


1 

y/2 


\F* 


z + 


Z + y/2 


< 


72 “ erf ( z > z - 


Now 


0? 


— dip 


1 — d I crf(z)z + 




< 1 — d(f). 


+ ^erf(z) 2 


= ^-^(1 - dtp) + ^erf(z) 2 > 0. 

y 7T Z 

From the above, 0 < h < 1 . 

The key observation is that for any fixed K and 2, F attains its maximum at c = 0 or c — 1. To show 
this, we explicitly compute the second derivative of F with respect to c, 

C, _ d ' 2 ^ 2 {K + 1){1 - dtp)h' 2 d ,2 crf(z) 2 

T 1 ‘1 I 


1 — d(j) 


2 h 2 


2h{\ -d^) 3 ’ 


where ' denotes differentiation with respect to c. Since F" > 0 , F is maximum at c = 0 or c — 1 as claimed. 
Hence if F is negative at these values of c, then we can conclude that F is a negative function. 
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The first term in the Taylor series expansion of F for small c is 


F= -# + 0(c 2 ) 

and hence F is negative for all small c. We now restrict to the CFL limit c = 1. 

In Figure 4.1, we plot F for \z\ < 100 and 2 < A < 100. It is clear that F is a negative function. 

We now examine the asymptotic behaviour of F for large values of 2 . For c £ (0, 1) and large |z|, 

F-(l- C )ln(l-c) + ^y|(ln(l-c) + l) + 0^. 

Hence F < 0 for large \z\. 

In paper [11], the positivities for both p* > 0 and p’E* - > 0 under the CFL condition have been 

proved. So, a distribution function f* with /* > 0 for the state {p*,m*,E*) can be constructed. 

Next we show that entropy increases in the second step where gas moves into cell Cj from its neighboring 
cells. It is sufficient to show only the case when gas from the left cell moves into cell Cj. Denote the 
quantities after one time step by {pj.m^Ej). See (3.7). 

Before proving the entropy condition in the above process, from Jensen’s inequality, it can be shown 
that pj > 0 and pjEj ifi- — 0, which means that the state ( pj,rhj,Ej ) satisfies the positivity condition. 
So, a gas distribution function f with / > 0 can also be obtained from this state. 

Lemma 3.2. Assume that pj,fhj,Ej are computed by (3.7). With the CFL condition, the entropy 
condition is satisfied in the process to obtain (pj,rhj, Ej). 

Proof. After some algebra, we have 

/ 9j In 9jdud( -crl ug r x In g 3 - X dudfi = \p. x F, 

J-oo Ju> o 2 |-z| + V7/2 

where 


F = <p 


In 


C(j) 

2(RT7t72), 


+ (K+ l)ln 


yfp 


1 = 

+ 


+ 


K + l . 


2v / tF ’ 


6 = z crfc(— z) + 


fix 


-2z 2 


i e “ , r / \ e z erfc(-2 ) 2 

= h 2 erfc(— z ) — — — . 

7T y/ix 2 


The goal is to show that the entropy condition is satisfied or equivalently, F(z, K ) c) < 0. As previously, 
we have introduced the non-dimensional number z = y/\ ~[Uj^ x and the number c is as defined in (3.13) 
but with j changed to j — 1. First note that O > 0 for all real values of z. To show this, note that 
dfi/dz = erfc(-z) > 0 and thus the minimum of 0 occurs at z = -oo where <j> = 0. Hence it is apparent that 
among the possible values of c € (0, 1], F is maximum at c = 1. Thus it is sufficient to demonstrate that 
F < 0 for c = 1. We shall assume this value of c for the remainder of this proof so that F is now a function 
of K and z. 

It can be shown that ip is negative for all z. Now for a fixed z, the term 


In 




2(| z | + v / t 72) 
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in F is maximum when 7 = 1 or K = 00 . The second term 


(K + 1 ) In 



is a decreasing function of K. This can be shown by taking its derivative with respect to K and it is 

D 

where 


■- ln(l + j/) + - — — — , 
2 V yJ 21+2/’ 


(K+ 1)0 2 ' 


Note that — 1 < y < 0. The derivative D can be shown to be negative for all y G (—1,0). Thus the second 
term achieves its maximum at K = 2 . Hence we conclude that 


F < <t> 


In 


0 


2(\z\ + 


- - In 
2 


( 1 + 5 $ 0 . 



G(z). 


For z G (0. 00), G z < 0 and since G(0) = —.5775 • * •, we have shown that G < 0 on [0, 00). For z < 0, G is 
maximum at 2 = — 00. As 2 — ► — 00, the first term of the asymptotic expansion of G is 


3e In |z| 
2y/n 2 2 


and so it is a negative function for 2 < 0. Thus we conclude that F is negative and thus the entropy condition 
is satisfied. We have finished the proof of the lemma. 

We plot F(z, K) in Figure 4.2. 

As a result, we have 


(3.14) 


Similarly, we have 
(3.15) 


/ gj In gj dud£ < a / ugj _ 1 In g 3 _ 1 dud^ . 

-00 Ju > 0 


/ gj In gjdud£ <—& ug \n gj + \dud^. 

-00 Ju < 0 


for the particles coming from the cell j 4 - 1 on the right hand side. 

After all terms of (p*, m*, £**), (p, (p, m, F 1 ) arc obtained, the flow variables in each cell Cj 

are updated according to Eq.(3.9). Since positivity is satisfied for each species (p*, m*, E*), (p, m, E) and 
(p, m,E), the distribution functions g*,g,g satisfy the conditions g* > 0,p > 0,p > 0 . In the collisional 
step, different species with its individual identification W* 9 W and W are mixed to form equilibrium states 
g*',g f and g f with a common velocity U and temperature A. In the collisional process, the individual mass, 
total momentum and energy are conserved, and the individual equilibrium states become 

K±l 

e _ A (( tt _c/) 2 +e a ) j 



(3.16) 


9 1 = P 



K + l 

e -\(( U -u) 2 +e) , 


9 



9 =P 


-\(( u -uf+e) 


where A and U arc determined from the total momentum and energy conservations Eq.(3.9), 

(p* 4- p + p)U = m* + m 4- m 

and 

(p* + P + p)(\u 2 + =E* + E + E. 

Lemma 3.3. The collision stage from (g*,g,g) to {g* f , g f ■> g') satisfies the entropy condition. 
Proof. Since 

9*> 0 , g> 0 , g> 0, 

and the individual mass, total momentum and energy conservations are satisfied, we have 

J 3 ^ n 9 dud^ 4- J g f In gdudf 4- J g Ingdud^ — J g* In g*dud£ ~ J g In gdud^ — J gin gdud£ 
~ J (9 ~~ 9 )l n 5 dudf J g ln(p /g*)dud£- f- J (g f — g) In gdud^ + J g\n(g f /g)dud£ 

+ J (9 ~ 9) In gdud£ + J g ln(g' / g)dud£ 

= J 9 ln(g* / g*)dud£ + J g\n(g / g)dud£ 4- j g ln(g'/g)dud£ 

— J 9 ( 9 / 9* ~ + Jg{g/g — l)dud£ 4- J g{g r /g — l)dud£ 

= J (g* - g*)dudf + j (g ~ g)dudE, + J (g f - g)dud£ 


In conclusion, we have 


(3.17), 


9 9 dud£ + j g'\ngdud£ + J g'kigdud£< J g* In g f dud£ -f J glngdud£ + J g In gdud£. 

Once wc have g ,g and g' . the total entropy of the distinguishable particle system inside cell Cj is 

H' — J g"' \ng*'dud£ + J g'\ng'dud£ + j g’ lng’dud£. 


and the total distribution function is 


9 = 9*' + 9 > +9' 


(3.19) = 




,++ 'L 


-\au-u) 3 +e) 


= ( P * + p + p) 


-H(u-u) 2 +s 3 ) 


With the updated (p, fn. E) inside cell Cj in Eq.(3.9), the total entropy If- 1 1 is composed of the sum of the 
individual entropies of three species, 


H n+1 = H' 


J g*' In g*'dud£ + J g' In g'dud^ + J g'Xng'dud^ 


(3.20) = P* In p* + p* (In — - l) + plnp + p^^(ln- -l)+plnp + p^^(ln- - 1). 

£ 7T Z 7T 2 7T 
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With the Lemma(3. 1-3.3) and the total entropy of three species at step n + 1, we have 
THEOREM 3.1. The entropy condition (2.4) is satisfied in the KFVS scheme. 

Proof. From Equations (3.11), (3.14), (3.15), (3.17), and (3.19), the new total entropy for the three 
species at cell j is 


H n+1 = R t 


= J g* f \ng* f dud£ + J g In gdud^F J g'lng'dud^ 

< f „• ^dud( + / 9 + f S (Lcm». 3.3) 

<HJ+ ^ (G? _ 1/2 - <^ 1/2 ). (Add Eqns.(3.11), (3.14) and (3.15)) 


Remark: the flow variables inside cell j at n + 1 do consist of three distinguishable species. 

For any numerical scheme, basically we are only remembering the conservative quantities inside each cell 
and the entropy is a function of the conservative variables. However, beside this, the entropy concept is also 
related to the information. For example, the entropy is different for a gas composed of one single color and 
a gas composed of two different colors. Numerically, at the beginning of each time step, we divide the gas 
into different cells. Consequently, the gases in different cells become distinguishable. For example, pf-i can 
be regarded as blue, pf as yellow and p™ +1 as red. As a result, inside cell C 3 at the end of time step n- hi, 
the gas p” +1 is composed of three species, i.c., red, yellow and blue, and the entropy Hj l+1 is the sum of the 
entropies of the individual species. The distinguishable effect of particles is purely due to numerical artifacts 
such as discretized space but they have a physical consequence. In order to remove the numerical effect at 
time step n + 1 inside cell Cj, we can numerically erase the different “colors” of the gas. More precisely, we 
can remove the individual history of the gas inside cell C 3 . As a result, the total density p CANNOT keep 
the information of the individual densities (p*,p, p), and the equilibrium state Eq.(3.19) goes to 


n+ 1 

j 


= P 



2 


e -X((u-U) 2 U 2 )' 


The corresponding entropy becomes 

H = J g In gdudf, 

(3.21) = plnp + p- ^ - h ln- _ i). 
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Note that for the same number of particles inside cell j at time n+ 1, there is quantitative differences in the 
entropies between distinguishable (Eq.(3.20)) and indistinguishable (Eq.(3.21)) system. This phenomena is 
related to the so-called Gibbs paradox [12]. 

The above post-process has no direct dynamical effect on the KFVS scheme in the updating of conserva- 
tive variables, and has no effect on the proof of the entropy condition in this paper. We are perfectly allowed 
to keep the individual species inside cell j at time n -F 1, and there is no need to take the above post-process 
to erase different colors and make them indistinguishable in terms of the updating conservative variables. 


4. Conclusion. The gas-kinetic scheme provides an approximate Riemann solution for the Euler equa- 
tions. The entropy condition for the Kinetic Flux Vector Splitting is proved in this paper. Based on the 
positivity and entropy analysis, we can conclude that the KFVS is one of the most robust schemes for CFD 
applications. 
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Fig. 4.1. Plot of F(z, K) at c = 1. 



Fig. 4.2. Plot of F(z,K). 
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